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ABSTRACT 

Hie  principal  assumptions  U6ed  by  R.  Kraichnan  in 
formulating  a  theory  of  turbulence  are  reviewed,  and  their 
meaning  is  examined  in  the  light  of  a  mathematical  model. 

A  regression  function  is  defined  and  used  to  evaluate  the 
triple  correlations.  The  importance  of  the  fourth  order 
cumulants  is  discussed.  Hie  merits  of  the  theory  are 
illustrated  by  comparison  with  experimental  results— in 

the  case  of  the  mathematical  "odel  —  of  grid  turbulence 

* 

(skewness  factor)  and  of  turbulence  at  very  large  Reynolds 


number. 
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I .  INTRODUCTION 


The  object  of  this  review  Is  to  examine  the  general  method  proposed  by 
R.  Kralchnan  several  years  ago.  Ve  will  probe  the  nature  of  the  assumptions. 

The  numerous  mathematical  details  will  receive  little  attention.  The  theory 
was  first  proposed  in  the  context  of  MHD  turbulence,  where  comparison  with 
experiment  was  not  possible.  The  problem  was  then  reduced  to  the  case  of 
ordinary  turbulence  (homogeneous.  Isotropic,  and  Incompressible),  and  some 
agreement  was  found. 

There  Is  little  doubt  that  Kralchnan' s  general  approach  can  be  extended  to 
MHD  problems,  to  turbulent  shear  flows,  or  to  plasma  turbulence.  In  fact,  a 
whole  category  of  problems  of  physics  may  benefit  from  any  success  In  the  field 
of  nonlinear  random  problems.  For  example  the  method  may  apply  to  nuclear 
magnetic  resonance. 

The  method  of  Kralchnan  Is  characterized  by  the  use  of  a  regression  function 
that  describes  the  average  effect  of  a  small  perturbation  imposed  on  the  tur¬ 
bulent  flow.  In  order  to  Introduce  this  function  as  clearly  as  possible  and  to 
show  that  it  can  be  measured,  I  shall  refer  to  a  model  problem.  This  problem 
Is  only  a  mathematical  exercise,  void  of  any  physical  reality,  but  easier  to 
treat  than  the  simplest  case  of  turbulence.  The  model  problem  has  the  advantage 
that  It  can  be  Integrated  numerically  with  a  large  computer  so  that  the  theory 
can  be  compared  to  experimental  results. 


II.  STATEMENT  OF  THE  PROBLEM 


Let  us  start  from  the  equations  of  Navier-Stokes  written  in  the  following, 
somewhat  Informal  form: 


(dv#/dt)  +  Z  v0(dv0/dx0)  +  (dp/dx.)  *  v(d2v#/&x#dx#) 
Z  (3v#/ax#)  »  0- 


(1) 


where  v(x,t)  is  the  Euler lan  velocity  and  p  the  pressure.  The  density  is  unity, 
v  is  the  viscosity,  and  the  dots  represent  assorted  indices  1,  2,  3»  The 
summation  extends  over  three  terms,  and  there  is  a  total  of  four  equations. 

In  order  to  obtain  a  system  of  ordinary  nonlinear  differential  equations, 
instead  of  partial  differential  equations,  we  Shall  now  use  a  Fourier  trans¬ 
formation  in  the  space  coordinates.  The  Fourier  transforms  u(k,t)  are  defined 
such  that 


v  *  Z  u#(k,t)  exp(ik#x,  )  (2) 

The  symbol  k  represents  a  wave  vector,  and  the  boundary  conditions 
correspond  to  a  periodic  flow.  The  interesting  part  of  the  flow  is  confined 
to  a  cubic  box  of  size  L;  therefore  the  lowest  magnitude  of  the  wave  vector 
shall  be  2n/L.  As  we  consider  higher  modes,  the  magnitude  of  k  increases  in 
Bteps,  and  many  different  orientations  are  possible.  In  principle,  the  magni¬ 
tude  of  k  could  go  to  infinity,  but  we  do  not  want  to  deal  with  wave  lengths 
comparable  to  microscopic  scales,  such  as  the  mean  free  path.  Thus  we  shall 
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specify  that  the  magnitude  of  k  has  a  maximum.  A s  long  as  this  maximum  is  high 
enough,  say  ten  times  above  the  Kolmogoroff  threshold,  it  has  no  significance. 
The  summation  of  £q.  (2)  therefore  covers  a  large  but  finite  number  of  terms. 

By  entering  £q.  (2)  into  Eq.  (l)  and  collecting  the  terms  proportional  to 

4  Vv 

e  ,  we  obtain  the  equations  of  Kavier-Stokes  in  the  following  form: 

dv^lO/at  -  Z  a  1Jh  Uj(k')  uh(k*;  -  vk2^  (3) 

The  coefficients  a^  are  exactly  defined  from  Jae  Havier-Stokeo  equations. 
The  two  Fourier  components  multiplied  by  a^  correspond  to  two  wavenumbers  k* 
and  k”  such  that  k*  +  k"  *  k.  Thus  the  summation  is  extended  over  all  triangles 
in  a  wavenumber  space.  Perhaps  one  could  formally  inelvde  all  pairs  of  vectors 
and  decide  that  a^  is  zero  unless  the  corresponding  vectors  form  a  closed 
triangle. 

For  each  magnitude  and  orientation  of  k  there  are  three  such  equations,  so 
that  we  have  now  a  very  large  system  of  nonlinear  ordinary  differential 
equations.  The  continuity  equation  is  implicit. 

I  will  now  recast  this  system  in  a  form  Which  is  more  appropriate  for  a 
review.  The  new  form  will  also  be  more  general  and  will  deal  only  with  real 
quantities.  Indeed,  every  u^(k,t)  and  every  a^  is  complex.  Beginning  vrth 
the  lowest  magnitude  of  k  and  a  first  orientation,  we  define 


u^k^t)  -  F1(t)  +  iFg(t) 

UgOtyt)  -  F3(t)  +  lF4(t)  (4) 

u^k^t)  «  F5(t)  +  iFg(t) 

Proceeding  to  the  next  orientation  or  magnitude  of  k,  we  write 

UjCk^t)  =  F?(t)  +  iFg(t) 

UgC^t)  -  F9(t)  +  iF10(t)  (5) 

UgO^t)  -  Fu(t)  +  IFjgCt) 

and  so  on  down  the  list  of  wavenumbers.  The  last  expression  is 

u_(k  ,t)  •  .......  +  1FW  (6) 

3  max  ■  R  '  7 

This  defines  the  real  quantities  F^  to  F^  where  H  is  a  very  large  number. 
For  the  simplest  case  of  turbulence  R  would  be  of  the  order  of  10^.  Each  F^ 
varies  with  the  time.  Rote  that  there  are  six  F’s  for  each  wavenumber,  so 
that  the  large  wavenumbers  correspond  to  functions  of  large  indices.. 

The  equations  of  Ravier  and  Stokes  can  be  written  in  terms  of  the  new 
functions  as  follows: 

ayat  - 1  =1Jk>yk  *  £  Vj  (7) 
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where  the  coefficients  axe  purely  real  and  related  to  the  coefficients 

a  .  of  the  previous  formulation.  The  coefficients  b  .  come  from  the  viscous 

1  jK  1 J 

terms.  These  equations  are  of  a  very  general  form,  and  they  encompass  a  wide 
assortment  of  problems.  Shear  flow  turbulence,  MHD  turbulence,  or  plasma 
turbulence,  which  do  not  contain  nonlinear  terms  worse  than  bilinear,  can  be 
reduced  by  Fourier  transformation  to  a  system  of  equations  of  the  form  of 
Eq.  (7). 

If  the  terms  in  c  .  are  negligible,  we  have  a  classic  linear  problem.  For 
ijK 

small  values  of  the  nonlinear  term,  a  method  of  perturbation  can  probably  render 
some  services.  In  turbulence,  the  nonlinear  terms  become  dominant,  and  a  funda¬ 
mentally  different  situation  appears.  In  fact,  the  terms  in  b^  either  become 
negligible,  or  add  only  some  minor  complications  to  the  general  strategy.  For 
simplicity,  we  shall  now  drop  these  linear  terms.  They  could  be  carried  along 
in  a  detailed  treatment. 

The  coefficients  c  grow  with  the  first  indices.  Indeed,  the  basic 
ijk 

equations  have  terms  such  as  vo(dv ^dx^)  that  lead  to  coefficients  proportional 
to  the  wavenumbers.  Thus  the  coefficients  having  large  first  indices  will  have 
large  values,  since  they  correspond  to  large  wavenumbers. 

It  follows  that  a  function  such  as  fluctuates  faster  than  a  function 

such  as  F^g  >  this  introduces  important  statistical  differences  between  the 
various  functions. 

The  problem  is  to  determine  the  statistical  properties  of  the  F^,  starting 
from  the  set  of  coefficients  generated  by  the  equations  of  Navier-Stokes, 

and  from  some  suitable  assortment  of  initial  conditions. 
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III.  A  MATHEMATICAL  MODEL 


In  order  to  examine  the  reasoning  of  Kraichnan  with  maximum  clarity,  I 

shall  occasionally  refer  to  a  model  problem,  which  is  not  related  to  any 

particular  flow  but  is  simply  defined  by  Eq.  (7)  with  a  convenient  choice  of 

1 

the  coefficients  c^^.  For  the  model,  these  coefficients  are  simply  taken 
from  a  collection  of  random  numbers  having  a  root  mean  square  of  unity  and  an 
average  of  zero.  The  initial  conditions,  that  is,  the  values  of  all  the  at 
t  =  0,  are  also  chosen  at  random. 

If  the  number  of  functions  F  is  less  than  100  or  so,  the  problem  can  be 
integrated  numerically  with  sufficient  accuracy  at  an  acceptable  cost.  I 
have  found  that  statistical  reproducibility  occurs  already  when  N  is  larger 
than  20,  and  satisfactory  numerical  experiments  can  be  carried  out  with  N  =  48. 
Ihis  is  a  great  relief  from  turbulent  cases  and  it  is  principally  due  to  the 
fact  that  all  functions  F^  have  similar  behavior  because  all  coefficients  have 
the  same  order  of  magnitude. 

The  state  of  the  flow,  or  that  of  the  model,  is  specified  at  any  time  by 
the  N  functions  F^.  Hiis  corresponds  to  a  point  in  a  space  of  N  dimensions. 

The  kinetic  energy  of  the  flow  is  defined  as  E  =  ^  F^F^,  and  it  follows  from 
Eq.  (7)  that  E  is  a  constant,  since  the  viscous  terms  have  been  discarded.  In 
the  case  of  the  model,  it  is  easy  to  adjust  the  coefficients  c^^  so  that  E  is 
also  an  invariant.  Then,  as  the  time  increases,  the  representative  point  stays 
on  the  surface  of  the  hypersphere,  in  the  N-dimensional  space. 

Starting  from  some  initial  conditions  at  t  =  0,  we  can  imagine  the  tra¬ 
jectory,  perhaps  like  the  solid  line  shown  in  Fig.  1. 
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Figure  1.  Solution  of  a  Turbulent  Energy- 
Preserving  Problem  in  a  Multi¬ 
dimensional  Space 


Figure  2.  Behavior  of  a  Component  in  a  Numerical 

Experiment 


IV.  THE  AUTOCORRELATION  FUNCTION 


The  physicist  confronted  with  the  signal  of  Fig.  2  immediately  realizes 
the  futility  of  a  theory  giving  every  wiggle  of  the  solution.  Some  statistical 
reduction  is  necessary,  and  the  first  question  is,  "What  is  the  power  spectrum?" 
The  fluid  dynamidst  will  ask  the  equivalent  question,  "What  is  the  auto¬ 
correlation?" 

This  function  is  defined  as 


H(t)  -  F(t)  F(t  +  t)  (8) 

In  the  case  of  turbulence,  this  function  varies  with  the  index  1  and  with 

the  corresponding  wavenusfcer.  This  function  H(k,T)  could  be  measured,  but  this 

2 

has  not  been  done.  Perhaps  the  results  of  Pavre  et  al.  could  furnish 
such  information. 

In  the  case  of  the  motel,  H  is  Independent  of  the  index  and  the  computer  can 
easily  determine  H(t).  The  results  give  the  bell-shaped  curve  shown  in  Fig.  3. 
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V.  EQUATION  FOR  THE  AUTOCORRELATION  FUNCTION 

The  theoretician  is  now  challenged  to  present  an  expression  for  H(t).  His 
first  step  will  be  to  form  an  equation  for  H.  For  this  purpose  we  denote 
F(t*)  as  F* .  After  multiplying  the  basic  Eq.  (7)  by  F'  and  averaging,  we 
obtain 

-(dH/dr).  ^rjfdFj/dt)  -  £  c1Jk  Fj'FjF*  (9) 

On  the  right  hand  side,  we  have  a  large  number  of  terms,  each  one  a  triple 

correlation.  Since  the  left  hand  side  is  of  the  order  of  magnitude  of  unity, 

in  some  appropriate  units  of  time,  each  triple  correlation  must  be  fairly  small 

Thus  the  various  components  F1  are  almost,  but  not  quite,  independent.  The 

classical  method  would  now  lead  to  equations  for  the  triple  correlations,  in 

terms  of  quadruple  correlations,  followed  by  some  drastic  assumptions  on  the 

3 

quadruple  correlations.  Kralchnan  attempts  to  directly  evaluate  the  triple 
correlations  by  means  of  a  regression  function.  We  shall  now  follow  this  new 
approach. 


Figure  U.  Effect  of  a  Small  Perturbation 
on  a  Turbulent  Solution 


Figure  5*  Effect  of  a  Small  Perturbation  on  a  Component  of  the 

Numerical  Model 


VI.  THE  REGRESSION  FUNCTION 


Let  us  come  back  to  the  model  and  examine  the  stability  of  the  solution. 
After  obtaining  the  solution  for  a  given  set  of  coefficients  and  a  given  set 
of  initial  conditions,  let  us  reset  the  computer  and  repeat  the  numerical 
integration.  Since  every  element  of  the  computer  is  fully  controlled,  the 
integration  is  completely  reproducible.  Even  the  errors  (truncation,  finite 
step  size,  etc.)  occur  in  exactly  the  same  fashion. 

At  time  t^  ve  can  stop  the  computer  and  modify  one  of  the  48  functions. 

Let  us  say  that  we  Increase  by  c  *  0.01.  (The  root  mean  square  of  is 
unity).  The  computer  is  then  instructed  to  resume  integration.  The  effect 
of  the  disturbance  is  illustrated  in  Fig.  4.  At  time  t^,  the  trajectory  makes 
a  sharp  step  in  the  direction  of  the  fifth  axis;  thereafter  it  follows  the 
undisturbed  trajectory  fora  while,  gradually  drifting  away. 

After  a  long  time,  the  perturbed  solution  is  far  away  from  the  undisturbed 
solution,  somewhere  on  the  hypersphere.  All  correlation  between  the  two 
trajectories  has  been  lost.  The  two  functions  F.  and  F,  can  be  plotted,  as 
in  Fig.  5.  The  difference  Af  ■  FK  -  Fc  can  be  divided  by  «,  so  that  it  Jumps 
from  zero  to  unity  at  t  ■  t^. 

After  a  long  time,  this  normalized  difference  becomes  very  large,  of  the 
order  of  l/«,  but  the  probability  that  it  will  be  positive  is  equal  to  the 
probability  that  it  will  be  negative. 


Once  again,  the  situation  calls  for  a  statistical  reduction.  Let  us 
repeat  the  procedure,  changing  only  the  Initial  conditions  at  t  ■  0.  For  each 
set  of  Initial  conditions,  we  have  another  normalized  difference,  Jumping  from 
zero  to  unity  at  and  eventually  becoming  very  large. 

A  family  of  such  results  Is  shown  In  Fig.  6.  Clearly,  the  behavior  of 
the  difference  can  be  predicted  for  some  time  after  t^,  but  the  long-range 
performance  is  completely  uncertain. 

The  computer  can  be  programmed  to  determine  the  average  of  the  family  of 
curves  shown  in  Fig.  6.  This  average  will  be  noted  as  G(t)  with  t  *  t  -  t*. 

For  K0,  G  ■  0  (see  Fig.  7).  This  is  the  regression  function,  and  it 
indicates  the  probable  effect  of  a  small  pulse  applied  at  t  «  0.  The  energy 
of  the  pulse  Is  gradually  distributed  among  all  the  N  functions  F.  Thus,  G 
describes  not  a  dissipation  of  energy,  but  a  loss  of  information. 

At  t  ■  tj,  the  probability  distribution  of  F^  and  of  any  other  component  is 
Gaussian.  Inmed lately  after  the  perturbation,  the  probability  distribution  of 
F,.  is  a  Gaussian,  shifted  by  the  amount  c.  In  time,  this  probability  will 
return  to  normal;  the  function  G  describes  this  regression  process. 

In  order  to  be  sure  that  the  individual  differences  indeed  become  very 

4m 

large.  In  absolute  value,  I  have  asked  the  computer  to  determine  the  function  K 


defined  as 


Sane  results  are  shown  in  Tig.  8,  for  various  values  of  IT.  The  horizontal 
scales  are  suitably  changed  so  that  the  functions  G  all  look  alike.  The 
functions  K  Indeed  grow  and  eventually  level  off  at  1 /«.  When  G  Is  almost  zero, 
the  growth  of  K  Is  exponential. 

As  IT  Increases,  note  that  K  follows  G  for  a  longer  tine  before  taking  off. 
Ohis  means  that,  as  the  system  becomes  more  intricate,  the  actual  effect  of  a  • 
pulse  can  be  predicted  vith  Increasing  accuracy.  Indeed,  the  gap  between  K 

and  G  is  a  measure  of  the  statistical  variations  between  differences  AF. 

23 

Far  IT  ■  10  ,  K  may  practically  vanish  before  the  final  ascension.  Such 
Systems  vould  have  a  "hydrodynamical"  response  to  perturbations! 
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VII.  APPLICATION  OF  THE  REGRESSION  FUNCTION 

Let  us  now  suppose  that  the  perturbation  Is  not  concentrated  on  one 
Instant  at  time  t^,  but  Is  continuously  applied  to  the  system. 

If  |(t)  Is  a  small  fluctuating  force,  generated  externally,  we  have 

(u) 

If  we  had  a  linear  system,  the  effect  of  the  external  force  could  be 
superposed  on  an  ordinary  solution,  leading  to  a  result  such  as 

t 

F  -SCt)  +  J  G(t-s)  f(s)  ds  (12) 

J  jo 

where  J  denotes  a  solution  In  the  absence  of  external  force.  For  a  linear 
system,  G  Is  simply  the  effect  of  a  single  pulse,  without  need  for  a  statistical 
averaging  or  a  limitation  to  small  amplitudes. 

With  a  nonlinear  system  of  sufficient  Intricacy,  we  can  perhaps  write  an 
equation  of  the  form  of  Eq.  12,  except  that  5  should  become  some  function 
closely  related  to  a  solution  of  the  undisturbed  equations. 

Essentially,  In  an  expression  such  as  Eq.  (12),  the  Integral  contains  those 
effects  of  t  which  can  be  predicted  by  means  of  the  regression  function,  while 
the  term  In  7  contains  the  effect  of  the  initial  conditions  plus  those  con¬ 
sequences  of  the  forcing  function  which  cannot  be  predicted,  ftius,  Eq.  (12) 
is  not  very  useful  If  we  look  for  a  formal  solution  of  the  problem. 
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However,  Eq.  (12)  becomes  valuable  when  we  form  the  correlation  t(t*)F^ 
between  the  forcing  function  and  the  solution  of  the  nonlinear  problem. 

From  Eq.  (12)  we  obtain 

TV  -  fV  +  /(Kt-s)  Y(t-s)  ds  (13) 

where  J  ■  ♦  't  Is  the  autocorrelation  of  the  forcing  function.  The  integral  is 
a  simple  operation  on  two  bell-shaped  functions.  The  correlation  la  far 
more  difficult  to  evaluate.  In  fact,  it  would  be  impossible  to  continue  the 
theory  if  we  could  not  dispose  of  this  termi  Essentially,  since  the  predictable 
effects  of  t  are  given  by  the  lntegral^telT^^lt  is  reasonable  to  assume  that 
♦ '  ?  is  negligible.  We  shall  do  so. 


vm.  THE  TRIPLE  CORRELATIONS 


The  previous  section  can  now  be  used  to  express  the  triple  correlations* 
Consider,  for  example,  the  equation  giving  FC;  and  assume  that,  among  the  many 

terms  of  the  right  hand  side,  we  find  one  proportional  to  the  product  F^Fg. 

i  i 1 

What  is  the  correlation  F^F^Fg  ?  By  pulling  one  term  outside  of  the  simulation 
and  indicating  its  presence  as  f  ■  c^g  F^Fg  we  have 


®5^dt  -  l  C5JkF3Fk  +  ♦ 


ONCTCHM  Lift 


The  argument  of  the  previous  section  can  be  used,  treating  one  term  as  a  small 
perturbation  applied  in  the  presence  erf  almost  all  the  other  interactions. 

It  leads  to  the  following  result: 


^gF^t)  -  c57g  J*  0(t-s)F^Fg  F7(s)  Fg(s)  ds 


This  approximation  is  valid  if  the  functions  F^Fg  and  F^  are  independent, 
except  for  the  linkage  by  the  coefficient  c^g.  Additional  terms  must  be 
entered,  proportional  to  Bimllar  effects  caused  by  and  Cg^7*  Such  terms 
are  labelled  "direct  interactions"  because  they  involve  only  the  three 
functions  forming  the  triple  correlation*  There  are  also  many  possibilities 
for  indirect  interactions,  involving  other  functions.  For  example,  a  term 
F^q  could  contribute  to  the  generation  of  F^  and  a  term  FgF^Q  could 


m 


contribute  to  the  generation  of  7300  •  Now,  If  the  product  ^qq^oo  <®Pear*  in 
the  equation  for  dF^/dt,  a in  approximate  analysis  shows  that  the  product 

affects  7^.  Thus,  through  the  action  of  7^  ve  have  another  linkage 
Kralchnan  assumes  that  the  direct  Interactions  are  the  dominant  ones.  This 


amounts  to  saying  that  the  various  Fourier  coqpoarats  are  correlated  three-by- 
three,  but  not  In  mare  numerous  groups. 


IX.  THE  EQUATION  FOR  H 


As  shown  In  Eq.  (15)*  a  triple  correlation  can  be  expressed  by  an  Integral 
over  a  function  0  and  a  quadruple  correlation.  Let  us  use  the  drastic  assumption 
of  MUllontschlkov  to  reduce  the  quadruple  correlation  to  a  product  of  two  double 
correlations.  Later  on,  In  Section  XII,  we  shall  comment  on  this  crucial  step 
and  examine  the  various  k*nda  of  quadruple  correlations  to  offer  some  Justifi¬ 
cation. 

Triple  correlations  can  now  be  given  in  terms  of  doubles.  Thus  the  road 
started  at  Eq.  (9)  leads  to  an  expression  of  the  type: 

T 

dH/dt  ■  C...C...J  Q(t— s)  H(s)  H(s)  ds  (l6) 


This  Is  the  first  of  two  master  equations  established  by  Kraichnan.  If  G 
Is  known,  Eq.  (l6)  gives  H. 


X.  THE  EQUATION  FOR  G 


The  application  of  a  perturbation  c  on  the  function  F,.  produces  the 
difference  AF,  seen  in  Fig.  6.  It  also  causes  a  disturbance  in  every  other 
component  of  the  motion,  and  ve  can  define  N-l  differences  6F,.  While  AF  is 
discontinuous  at  t  ■  t^,  the  quantities  6F^  are  continuous  but  each  shows  a 
discontinuity  in  some  derivative.  Thus,  if  F^  is  directly  modified  by  F,.,  the 
derivative  of  6Fq  will  jump.  If  F^  is  directly  affected  by  F^  but  not  by  F,., 
the  second  derivative  of  6F^  will  Jump,  etc. 

The  averaged  value  of  any  function  6F.  is  zero  (j  $  5)>  since  the  pertur- 
bations  are  "modulated"  by  various  functions  F^. 

Immediately  after  the  perturbation,  and  during  the  decay  of  G,  the  difference 
AF,.  and  the  many  other  differences  6Fj  are  small.  Thus  the  basic  equations  can 
be  linearized  and  we  write: 


d/lt  4F5  -  2  l  =5J16FjFk 


(17) 


d/dt  6F.  -  V  c,K  AFCF  +  2  Y  c.6F  F 
^  J  L  J5P  5  P  Li  Jpq  p  <3 


(18) 


In  order  to  determine  G  *  AF^/c,  we  must  eliminate  the  quantities  6Fj.  The 

equations  for  6F  show  that  every  one  is  a  response  of  the  system  to  a  per- 
J 

turbation  AF,.  F^.  The  predictable  part  of  this  response  is 


XI.  THEORY  FOR  THE  MODEL 


In  the  case  of  the  mathematical  model,  the  second  master  equation  reduces 
to  a  form  Identical  with  the  first  equation,  so  that  one  finds  G  *  H,  provided 
that  T  is  positive,  Thus  the  single  equation  is 

dB/dT  -  -c...c...  f  H(T-s)  H^s)  ds  (21) 

0 

It  can  be  Integrated  easily  by  numerical  techniques,  and  the  results  are 
shown  In  Fig.  9»  The  solid  curve  gives  the  experimental  results,  the  dashes 
and  dots  show  the  solution  of  Eq.  21  from  the  theory  of  Kralchnan. 

The  dots  show  the  experimental  values  for  G.  The  discrepancies  between 
theory  and  experiments  are  probably  due  to  numerical  errors  In  the  integration 
of  Eq.  21  or  to  limited  averaging  samples,  or  to  the  fact  that  the  theory  applies 
for  N  ■  °°,  while  the  experiments  apply  to  N  ■  48. 

In  order  to  compare  the  new  theory  with  older  methods,  we  can  form  an 
equation  for  triple  correlations  in  terms  of  quadruple  correlations  and  use  the 
drastic  reduction  to  doubles.  The  result  is  shown  In  Fig.  9  by  the  thin  curve, 
which  plunges  to  -  09 .  Thus,  the  new  theory  marks  a  clear  progress. 


H  (THEORY) 


Heavy  solid  line:  From  numerical  experiments. 

Heavy  broken  line:  From  Kralchnan's  theory.  Ea  (21) 
min  solid  line:  From  discard  of  cumulantJ!  * 

OMn  broken  line:  Regression  function,  from  numerical  vorl 


Figure  9.  Autocorrelation  Functions 


XII.  HOLE  OF  THE  QUADRUPLE  CUMULANTS 


In  order  to  examine  in  detail  the  difference  between  old  and  new  theories, 
we  shall  start  by  establishing  an  exact  relation,  a  refreshing  episode. 
Multiplying  the  basic  Eq.  (7)  taken  at  time  t  by  the  same  equation  taken  at 
time  t*  and  averaging  leads  to 


(aydtKd^/dt)  -  ScyjE clpqFJFJPpPq  (22) 

If  the  process  is  statistically  stationary  in  time,  the  first  integration 
gives 

(23) 

0 

Note  the  double  summation.  The  right  hand  side  of  Eq.  (23)  can  be  regarded 

as  a  product  between  a  matrix  c,  ..  c .  and  a  matrix  formed  by  various  quadruple 

ljK  lpq 

correlations . 

Let  us  now  approximate  these  quadruple  correlations  by  products  of  doubles. 
Since  FjFj  vanishes  unless  i  ■  J,  the  only  terms  of  Eq.  (23)  that  do  not 
disappear  are  those  found  along  certain  diagonals  of  a  matrix. 

Thus  the  exact  Eq.  (23)  reduces  to  an  expression  typical  of  these  theories 

dH/dT  »  -  )  c...c...  J*  HH  ds  (24) 

^  0 
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1 

1 


1 


Note  that  the  summation  has  been  reduced  to  a  single  dujj:.  If  there  is 
a  small  error  in  the  evaluation  of  each  quadruple  correlation  along  the 
diagonals,  it  will  have  a  small  effect  on  the  result.  However,  if  there  is 
a  small  error  on  each  quadruple  correlation  away  from  the  diagonals,  the  very 
large  multitude  of  these  terms  may  produce  significant  errors.  Urns,  it  is  the 
neglect  of  the  cumulants  off  the  diagonals  that  is  likely  to  cause  trouble. 

The  theory  of  Kraichnan  leads  to  Eq.  (l6),  which  reduces  to  Eq.  (24)  if  we 
take  G  •  1.  Thus,  the  presence  of  the  regression  function  in  the  final  results 
is  significant.  Since  G  represents  some  effects  of  the  complicated  coupling 
between  the  modes,  it  perhaps  accounts  for  the  off-diagonal  terms.  How  well 


is  not  known. 

Finally,  if  we  return  to  Eq.  (l6),  we  remember  that  Kraichnan  uses  the 
simplified  form  for  the  quadruple  correlations  at  some  stage  of  the  argument. 
However,  this  reduction  involves  only  special  correlations  along  the  diagonals 
and  does  not  abridge  a  double  sum  to  a  single  sum.  Thus,  it  seems  that  the 
cumulants  can  be  neglected  only  along  the  diagonals. 


XIII.  APPLICATION  TO  GRID  TURBULENCE 


The  turbulence  produced  by  a  grid  placed  in  a  wind  tunnel  is  generally 
characterized  by  a  low  Reynolds  number.  With  a  mesh  Mq  and  a  free  stream  UQ* 

we  shall  refer  either  to  Rj^  ■  Uc of  the  order  of  15*000,  or  to 

4-  **■ 

R ytv  (O.IR^)  *  varying  between  20  and  70. 

The  energy  spectrum  is  too  narrow  to  display  the  Kolmogoroff  range. 

Kraichnan  applied  his  theory  in  the  case  R^  ■  M),  assuming  a  convenient  form 

for  the  shape  of  the  spectrum  at  t  «  0.  The  turbulence  decays*  so  that  the 

functions  H  and  G  vary  slowly  with  the  time.  The  successive  spectra  are  shown 

in  Fig.  10.  Note  the  decay  of  the  energy  and  the  production  of  small  vortices. 

The  decay  proceeds  according  to  a  similarity  law. 

For  this  flow,  Kraichnan  determined  the  skewness  factor  S  defined  as 

S  -  (dv^/dx^3/  (25) 

where  v^  and  x^  are  along  the  mean  flow  direction.  Note  that  S  is  nondlmenslonal 

This  factor  plays  an  essential  role  in  the  study  of  nonlinear  turbulent  processes 

It  is  related  to  short  range  triple  correlations  and  it  is  also  proportional  to 

the  rate  of  production  of  the  mean  square  vcrticity,  and  to  a  rate-of-deform- 

5 

atlon  parameter.  Certain  theories  are  based  on  an  expression  for  the 

6 

eddy  viscosity*  which  introduces  one  nondlmenslonal  universal  constant 


that  must  be  matched  with  experimental  data.  Oils  constant  is  none  but  the 
factor  S,  in  a  special  form. 

In  Kraichnan's  treatment  It  is  necessary,  at  t  *  0,  to  specify  the  triple 

correlations;  he  started  from  provisions  for  zero  triple  correlations,  which 

corresponds  to  S  ■  0  in  Fig.  11.  The  rapid  initial  growth  of  S  is  compatible 

7 

with  previous  work  of  Proudman  and  Reid,  which  used  the  drastic  treatment  of 

quadruple  correlations.  However,  the  new  theory  finds  that  S  rapidly  levels 

off  and  retains  the  value  S  ■  0.4  during  the  decay.  This  is  in  excellent 

agreement  with  the  experiments  of  Townsend  and  others. 

Note  that  Kralchnan  does  not  adjust  any  constant  in  his  determination  of  S. 

He  also  treats  several  similar  cases,  always  with  the  limit  S  -*  0.4.  A  kine- 

5 

matlc  argument  shows  that  S  cannot  exceed  a  value  near  0.8. 
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Figure  11.  Skewness  Factor 

(According  to  Kraichnan  for  ■  40.  The  Value  Factor  0.4  Agrees 

with  Many  Experiments.) 


XIV.  THE  KOLMOGOROFF  SPECTRUM 


In  the  atmosphere,  the  ocean,  or  a  sufficiently  large  apparatus,  the 

spectrum  of  turbulence  shows  a  portion  falling  as  k-'^.  We  shall  consider 

the  case  of  stationary  turbulence,  so  that  we  must  Imagine  some  large  scale 

mechanism  furnishing  the  energy  necessary  to  maintain  the  flow,  despite  the 

4 

energy  dissipation  c. 

In  a  first  attenjpt  to  apply  his  theory,  Kraichnan  met  with  some  difficulties 
In  the  limit  of  very  large  Reynolds  numbers  (say  ■  300),  the  equation  for 
G(k,T)  takes  the  form  of  a  product 


dG(k,T)/dT 


|^A(k*,k")||  J^t(k,r,s)  dsj 


(26) 


where  the  Integration  in  s  contains  only  a  function  of  t,  k,  and  s.  The  sum, 
however,  contains  a  contribution  from  each  triangle  In  the  wavenumber  space 
such  that  k*  +  kw  *  k.  If  either  k*  or  k"  becomes  smaller  or  comparable  to 
the  wave  number  K  at  which  the  spectrum  Is  maximum,  the  behavior  of  A  requires 
special  attention.  As  a  result,  the  final  expression  for  the  energy  spectrum 
takes  the  form 


E(k)  -  c2/3  £(a/ k5/3)  +  (b/fe^k3/2]  (27) 


Figure  12.  Autocorrelation  Functions  Produced  by  a  Sinple  Shear  Wave 
Transported  at  Various  Uniform  Speeds 
(The  Averaged  Effect  is  Shown  by  the  Heavy  Line) 


The  term  proportional  to  a  comes  from  values  of  k'  and  kM  larger  than  K.  The 
term  proportional  to  b  comes  from  the  small  values  of  k*  or  k",  that  is,  from 
the  large  scale  vortices. 

It  is  clear  that  if  k  is  large,  the  above  expression  falls  as  k”^.  If  the 
singular  behavior  did  not  occur,  b  would  vanish  and  the  spectrum  would  fall  as 
expected. 

This  situation  prompted  Kralchnan  to  a  reexamination  of  the  first  form  of 
his  theory.  He  found  that  the  argument  that  I  have  outlined  in  the  previous 
pages  suffers  from  a  serious  defect.  It  can  be  corrected  if  the  correlations 
are  defined  in  a  sort  of  local  frame  of  reference  moving  with  the  local  fluid 
velocity.  The  papers  dealing  with  the  revised  theory  are  easily  identified: 
the  title  always  contains  the  name  of  Lagrange.  They  are  not  easily  read. 

To  give  you  an  idea  of  the  nature  of  the  difficulty,  I  shall  consider  a 
situation  so  simple  that  there  are  only  two  modes.  First  let  us  imagine  a 
fluid  at  rest,  animated  by  a  single  stationary  shear  wave  such  that 
v(x,y,z,t)  •  sin  ax. 

Clearly,  the  correlation  v(x,t)  v(x,t  +  t)  is  independent  of  T  since  the 
shear  wave  is  stationary.  Let  us  now  displace  the  fluid  with  a  supplementary 
velocity  U,  which  is  constant  in  space  and  in  time.  We  assume  that  the  shear 
wave  is  transported  with  the  velocity  U. 

Then  the  velocity  observed  at  a  fixed  point  varies  in  time  and  the  time 
correlation  will  fluctuate.  For  a  particular  value  of  U,  the  correlation  will 
be  given  by  one  of  the  thin  lines  shown  in  Fig.  12.  If  we  change  the  magnitude 
of  the  velocity  U,  we  get  some  other  oscillating  correlation  function;  only  the 


period  varies. 

An  average  over  an  ensemble  of  different  velocities  U  or  over  time,  if  U 
changes  slowly,  will  yield  the  bell-shaped  correlation  function  shown  by  the 
thick  line  in  Fig.  12. 

Thus,  the  transport  of  a  single  shear  wave  by  an  assortment  of  large  scale 
motions  can  produce  a  certain  kind  of  correlation.  However,  one  should  not 
expect  that  the  theory  of  Kraichnan  will  succeed  in  this  case.  The  theory 
must  be  applied  only  to  those  correlations  that  are  not  produced  by  random 
translations  of  the  same  shear  wave. 

It  is  too  early  to  say  that  the  modification  proposed  by  Kraichnan  is  the 
only  possible  one,  but  it  is  not  surprising  that,  once  the  treatment  for  low 
k'  and  k"  ie  revised,  the  spectrum  falls  as  k-^^.  The  Kolmogoroff  law  is 
already  in  the  Eulerian  theory.  It  must  also  be  noted  that  the  Lagranglan 
theory  gives  a  numerical  value  of  the  universal  spectral  constant  a  that  is 
in  good  agreement  with  experiments. 

Turbulence  at  high  Reynolds  numbers  raises  many  interesting  questions. 

What  is  the  skewness?  It  has  not  yet  been  measured  experimentally.  One  could 
al6o  take  the  signal  proportional  to  a  velocity  fluctuation  and  pass  it  through 
a  simple  low-pass  filter.  The  skewness  of  the  filtered  velocity  could  be 
determined  and  perhaps  compared  with  theoretical  results. 


XV.  CRITICISM  OF  THE  THEORY  AND  COMffiNTS 

So  far,  the  theory  has  not  predicted  many  remarkable  properties  of  tur¬ 
bulent  flows,  except  in  the  field  of  MHD  turbulence,  where  the  confirmation 
may  take  a  long  time. 

In  ordinary  turbulence,  measurements  of  H(k,T)  would  be  interesting, 

since  a  comparison  would  be  possible.  A  special  effort  should  be  made,  perhaps, 

to  specify  the  odd  correlations  of  order  higher  than  three.  If  the  direct 

interaction  approximation  is  valid,  the  triple  correlations  should  determine 

those  of  order  five,  seven,  etc.  This  would  explain  the  results  of  Frenkiel 

8 

and  Klebanoff . 

Some  scientists  have  not  fully  accepted  the  cascade  theory  and  Richardson's 
poetry.  In  electronics-engineering,  the  energy  supplied  at  60  cps  is  often 
converted  to  much  higher  frequencies,  say  megacycles,  without  any  intermediate 
oscillations.  The  spectrum  has  essentially  two  separate  peaks.  If  the  high 
frequency  signal  is  intensely  modulated,  the  spectrum  could  fill  the  entire 
range.  Thus,  one  can  perhaps  conceive  other  theories  of  turbulence,  with 
direct  links  between  the  very  large  and  the  very  small  eddies.  Intexmittency 
might  perhaps  play  an  essential  role.  So  far  these  are  only  speculations, 
and  the  theory  of  Kraichnan  still  stands  as  the  best  that  we  have. 
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